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Abstract 

^-H We carry out a comprehensive analysis of models for top Aps at CDF in light of new top data arriving 

o 

CN from the LHC. We begin with a careful Tevatron analysis, considering in general which sets of effective 

^ vertices give rise to a large forward-backward asymmetry while suppressing the contribution to the total 
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tt cross-section. We show on general grounds that scalar models struggle to produce sufficient asymme- 
tries consistent with CDF observations, while vector models can produce a large asymmetry with a less 



r^ significant tension in the total cross-section and tt invariant mass distribution at the Tevatron. We examine 

the essential observables of these models for top physics at LHC7 with 1 fb ^ of data, including the total 



cross-section, invariant mass distribution and number of additional jets in tt events. In the case of t-channel 

mediators, the LHC total cross-section places a strong constraint on light mediators, while the Tevatron 
(N 
^ invariant mass distributions place strong constraints on heavy mediators that are able to produce the asym- 



metry. Heavy axigluons are becoming increasingly squeezed by LHC7 tt and dijet resonance searches. We 
conclude that LHC7 top analyses are rapidly closing the window for viable models of the CDF top ApB- 



Contents 
I. Introduction |2] 

II. Effective Vertices and Top Aps 13 

III. Comprehensive Search for Models at the Tevatron [15] 

IV. Implications for Top Physics at the LHC 
V. Conclusions 

A. Parton Level Asymmetries 

B. Event Generation 

1. Fast Detector Simulation and Object Reconstruction 

2. Event Generation and Parameter Scan Strategy 

3. Cluster Pipeline Setup 

References 

I. INTRODUCTION 

The Large Hadron Collider (LHC) is providing an unprecedented probe of top quark properties. 
While the Tevatron has to date collected on the order of a thousand tops, the LHC, with 1 fb^ of 
data has already nearly an order of magnitude more tops. The improvement is due both to a larger 
production cross-section, and to improved rapidity coverage for leptons in semi-leptonic and fully 
leptonic top analyses. In terms of the percentage error on the total cross-section, the 7 TeV LHC 
(LHC7) results are already competitive with the Tevatron [fTl with only 35 pb~^ [|2l, while the 
invariant mass distribution with just 200 pb^^ extends to a higher mu of 2.5 TeV [3] as compared 
to the Tevatron reach of 1 .8 TeV flU. At 1 fb^, the top quark properties will be far better measured 
than at the Tevatron. 

At the same time, the Tevatron as a pp machine is better able at the outset to measure a forward- 



backward asymmetry.^ The asymmetry in a particular invariant mass bin, mjf j, is defined by 

.tu . ^ NjAy > 0, mu,^) - NjAy < 0, m,,-,) 

^"""'^^ N{Ay > 0, m«-,,) + N{Ay < 0, m„-,) ' ^ ^ 

with Ay the rapidity difference between a top and an anti-top. The recent CDF anlaysis shows 
ApB = 0.475 ± 0.114 for ma > 450 GeV [7j| at the parton level (or ApB = 0.266 ± 0.062 at 
the signal level)^, while the Next-to-Leading Order (NLO) Standard Model (SM) predicts much 
lower values 0.088 ± 0.013 (or 0.043 ± 0.009 at the signal level) Ml^, corresponding to a 3Aa 
deviation (3.6a at signal level). A measurement of the asymmetry with fully leptonic tops has also 
been made which is roughly consistent with the measurement in the semileptonic channel [fT3l . 
The DO collaboration also observes a larger than predicted asymmetry [fT4l . 

Because this asymmetry is so large, any new physics (NP) that could generate such an asym- 
metry must have large couplings to the top as well as to the light quarks in the initial state. A very 
large number of models have been proposed in the literature. However, from a phenomenological 
point of view, these models mainly fall into only two categories according to the nature of the new 
particle exchange: (i) s-channel exchange of vector mediators (e.g. axigluon models) [fT5l - l35l or 
(ii) t-channel exchange of flavor- violating mediators [|36] - l52l ^. Comparative studies of these mod- 
els have also been carried out [|54l - [6T| , and their implications for top observables at the Tevatron 
observed [I6UI62II . The s-channel mediators often have maximally axial couplings (though there 
are exceptions such as [|35l ). while the t-channel mediators connect a light quark to the top quark 
in a way that appears to maximally violate flavor. A number of studies on the implications of these 
models for LHC physics have also been carried out [| 45ll63l -[66l. 

There are a large number of possibilities for the spin, color, flavor and electroweak representa- 
tion of a new field that fits into the two categories mentioned above. In the literature these have 
been mostly built and studied one by one. Here, by contrast, we are motivated to extract general 
features to determine which effective vertices are able to generate the large ApB while contributing 
a small amount to the total cross-section. We find that the form of the matrix element itself allows 
one to make general conclusions about which classes of models are successful in generating a 
significant asymmetry. 

' Though see 15] [6l for efforts to make a measurement of the forward-backward asymmetry at the LHC. 

^ Throughout this paper we use "signal level" to refer to background subtracted, raw measured quantities in the 

detector, and "parton level" to refer to unfolded results which attempt to subtract detector effects from the results. 
^ There is another class of models that can create effective axial QCD coupling from NP |53 1. 



We find on general grounds that perturbative^ scalar models typically can produce no more than 
a 10 — 20% "parton level" asymmetry for rritt > 450 GeV, which is only somewhat larger than 
the asymmetry produced in the SM (at ~ 9%) and well below CDF's parton level central value of 
~ 48%.^ The reason is simply the combination of the Mandelstam variables that enters into t- and 
u-channel processes for scalars; the statement is independent of the color (singlet, triplet, sextet 
or octet) or flavor representation of the state. By contrast, t-channel vectors have a matrix element 
that is conducive to producing a large asymmetry with a relatively small contribution to the total 
cross-section. 

We systematically enumerate the possibilities for the quantum numbers of t- and s- channel 
mediators that can produce an asymmetry and show that classes of models are strongly disfa- 
vored based on a small contribution to the total asymmetry or large contribution to the total ti 
cross-section at the Tevatron. This paper is intended to be a companion to our earlier paper on 
ApB ll6ni . which carried out a systematic comparison of NP models to the data. This was the 
only theory paper to carry out the full top reconstruction in order to compare results at the signal 
level. We found that there were large acceptance effects which changed the extracted parton level 
comparison between the SM and the NP models.^ 

With LHC data quickly arriving, however, the source of strong constraints is rapidly changing, 
and we are particularly compelled by the fact that the LHC collaborations are now analyzing 
unprecedented amounts of top data that will clearly rule out a large swath of models. We examine 
observables from the LHC, such as total cross-section, ti invariant mass distribution, and the 
number of additional jets in tt events. In order to carry out our analysis, we have done a systematic 
scan in mass and coupling space for a broad class of models, described in appendix |Bj For a subset 
of the models that give the best fit to the data we generate 5 million events, applying cuts and 
mif reconstruction mirroring the ATLAS analysis [3] in order to compare to LHC tt distributions. 
Many models will be strongly constrained by these analyses with just 1 fb^^.' In our earlier 

^ Scalar models with larger couplings can achieve larger asymmetries, though at the expense of a larger contribution 

to the total ii cross-section. 
^ The data-level asymmetry yields a result about a factor of two lower than the parton level result, which has been 

confirmed by the theoretical study of |61|. A comparison of a parton level theoretical result to the signal level 

asymmetry is not valid, and will underproduce by more than 2(7 the observed asymmetry. 
^ See also |50| for a discussion of the acceptance effect at the parton level. 
^ These statements must take into consideration, however, uncertainties in next-to-leading-order (NLO) corrections 

to the NP contributing to the total cross-section and invariant mass distribution. 



paper on searching for flavor-violating resonances at the LHC, we proposed top-jet resonances as 
a means to search for t-channel mediators [|63l . Such a search is complementary to the analysis 
here. Many t-channels models will be constrained by existing analyses, but the models that survive 
can have an imprint in top-jet resonances. 

The outline of this paper is as follows. In the next section we discuss the classes of models that 
could generate the forward-backward asymmetry at the parton level, examining the asymmetries 
that are generated by the possible effective vertices, and drawing conclusions about which classes 
of models are viable. The reader who is interested only in the numerical results can skip this 
section, and move on to Sec. |in| referring to Sec. ^ solely for a discussion of our conventions. 



In Sec. (Ill), we carry out a systematic scan of models at the Tevatron, choosing a set of models 



as benchmarks for simulation of the large data sets necessary for invariant mass distributions. In 



Sec. (IV), we then examine the expected top properties at the LHC for the classes of models we 
consider. In the appendices, parton level asymmetries, as well as a detailed discussion of our 
analysis pipeline, can be found. 

n. EFFECTIVE VERTICES AND TOP ApB 

Broadly speaking, either s-channel or t- (or M)-channel resonances can generate the top 
forward-backward asymmetries at tree level. We show the diagrams that contribute both to the 
Tevatron Afb and tt production at LHC in Fig. [T| The structure of the differential cross-section 
for models that produce the asymmetry through t-channel exchange of a top-flavor-carrying me- 
diator takes the same basic form according to whether the mediator is spin-0 or spin-1. Let the 
effective Lagrangian involving top and up quarks take the form 

_ r tigLPi + gRPR)KuM^ + H.c. spin-0 
I t^igLPi + gRPRKuM; + H.c. spin-1 

where t = t for singlets and octets and t = f^ for anti-triplets and sextets, and t° are the color gen- 
erators of a representation r — 3 x 3 Hermitian matrices that contain Clebsch-Gordon coefficients 
connecting two (anti-)quarks, normalized so that Tx(t'!;.t\) = ^5"'''. For singlets, we take t" = 1. 
Note that we are restricting ourselves to couplings to top and up quarks, though our results should 
not qualitatively change given couplings to down-type quarks. We also only consider single medi- 
ator production; double mediator production is only important for light colored states, which are 
not present for the models we consider. 






(a) s-channel qq (b) s-channel gg 

9 i . ^9 i 




q 



M 




(c) tt production 

9 M 




/ 



M 



q 



(d) t channel single mediator (e) s channel single mediator (f) u channel single mediator 

FIG. 1: Tree level production diagram involving the mediator M and the coupling gM- 

A large number of models of NP can generate these effective vertices, involving, in addition 
to the color group, potentially SU{2) and flavor representations. There are a limited number of 
flavor symmetric models that can generate the top ApB in the t channel while satisfying existing 
constraints. We show in Table |l] the possibilities. There are also interactions that connect Ql 
to Ql, but these models with flavor symmetries are typically highly constrained by light quark 
observables since they mix with SM CKM physics. We do not consider them further. Interactions 
connecting, e.g., Ql to ur through a spin-1 color triplet or sextet "diquark" are also possible, but 
as we will soon see (see Fig. [6]), the dominant t-channel interaction for top ApB does not give rise 
to a significant positive asymmetry. We refer the reader to [IFTl l68l for a complete tabulation in the 
scalar mediator case of the possible flavor symmetries and to Il46l |69l for discussion of Aps in the 
context of Minimal Flavor Violation. In any case, the general observations that we make on the 
basis of the effective vertices in Eq. (|2]) will be relatively independent of the flavor representation, 
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TABLE I: Flavor symmetric interactions (in schematic notation) involving at least one ur quark that can 
mediate a significant positive top forward-backward asymmetry in the i-channel. (See also ll69]| .) 



and we make the appropriate qualifications where necessary. For example, in flavor symmetric 
models, states in both the t channel and the s channel can contribute to the total asymmetry. 
Scalars in the s channel don't contribute to the forward-backward asymmetry, but can have an 
impact through their interference with t-channel scalars that do generate the asymmetry. 

Given the large number of possible combinations of s- and t-channel resonances from the flavor 
symmetric models, one despairs of ever being able to derive the characteristics of the state that can 
generate the asymmetry. However, we will find that in the t (or u) channel, the amplitudes have 
very distinctive shapes dependent on whether the state is a vector or scalar mediator particle. We 
examine these characteristic features, and use it to draw conclusions about the nature of the medi- 
ator from the invariant mass dependence of Aps- These conclusions are robust independent of the 
particular flavor symmetric model that one employs, and allows one to make general statements 
on the types of characteristics that are necessary for generating a large top ApB- 

The cross-sections arising from the NP interactions ([2]) and SM interactions are given by 

da /3 



dcos6 327rs 



(^SM + Ant + A 



sqj } 



where UMI EH [H 



A 



SM 



^1 

9 



l + c^ + -r^ 



ytji 






i9l + 9l) ^f" spin-1 



stM 



i9l + 9l) 



2\ tt+smi 



{i -H- u for diquarks) 



(3) 

(4) 

(5) 



spin-0 



Color rep: 18 3 6 

C(o) 4 -2/3 1 -1 
C(2) 9 2 3/4 3/2 

TABLE II: Color factors for color representations of flavor-changing mediators. 
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Aq = ^^| *M '' ^^^^'^ (t4->Mfordiquarks). (6) 

9 ^ (^i±^^ spin-O 
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Here C^Q^ and CLa are color factors depending on the color rep of the mediator.^ We have also 
defined 



ce = Pco^e (3 = ^l- Amj/s, (7) 

ti =t — mf Ui = u~ ml. (8) 

The Mandelstam variables are related to the scattering angle via 

i= -s{l-CB)/2 + m] and u = -s{l + ce)/2 + ml. (9) 

Note that we have not taken into account interference between NP contributions which can arise in 
flavor symmetric models. For example, s-channel flavor conserving and t-channel flavor changing 
diagrams may interfere. These new contributions do not give rise to any new types of terms 
(modulo mass terms in propagators) in the interference amplitude for the vector states, but do give 
rise to new contributions for the scalar states. We discuss these terms later, but suffice for now to 
comment that the new terms will not change our qualitative conclusions. 

A flavor-conserving vector can give rise to an asymmetry at tree level if couplings to top and 
up have nonzero axial parts. Given the NP interaction Lagrangian 

/:np = {qT^l^{glPL + glPR)q + t T^7^(^i^L + ^^P«)t) G'% (10) 

the scattering cross-sections calculated through these interactions are 1.56.1 : 

Ant = ^iT^^rS^n^ (^^^^' + ''* + 2"^*^) + ^"^^* - '^)) ' (ii> 

y b\bQ -t- iiiqL q) 



• Specifically, CJ„^ ^ -^Tr U^T'^t^T^] and CJ^^^ = Tr (i^4) Tr (i^4) where ^ = -1(1) and T^ = T'^{T^' ) 
for octets and singlets (anti-triplets and sextets). 

8 



and 

Aq = Ir^-^^^, (iaf + 91') (igi" + a'a^^l + i^) + 2g\g'^2mls) + g-g^{ul - il 

(12) 
where 

g^^{gl±gl){g\±g'^) and sc^s-ml. (13) 

For a color singlet rather than a color octet, the interference term vanishes and the squared term is 
scaled by a factor C(^2)/'^f2) = 9/2- 

We now assemble these results using the parton distribution functions to gain a strong quantita- 
tive understanding of which types of interactions can give rise to the observed forward-backward 
asymmetry. The cross-section for the process pp — )■ tt is given by: 

o{^s) = Tjij ds dx — / dcosO fi(x)fj j — J aij(cos9,s). (14) 

We define 

F,,is,s)= I dx-h{x)fJ—\ (15) 

Js/s -^ \SX / 

Then the differential cross-section as a function of parton energy s can be expressed as 

da(s) 1 

—— = -i:ijFij{s, s)d-ij{cos9, s). (16) 

dsd cos 9 s 

Of course, all the cos 9 dependence is in the parton-level differential cross-section. If only one 
kind of initial state parton contributes to the cross-section, then the PDF completely factors out of 
the differential forward-backward asymmetry, defined as a function of s by 

Afb{s) = ^ ^ ,. ^-+f-v (^^^ 



where 
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/ dz{crij{z,s) ±aij{-z,s)) . (18) 
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Suppose we are interested in a NP model with a nonzero contribution to the cross-section term 
generated through uu — )> ti. We may then write the forward-backward asymmetry as a function of 

\/S as 
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FIG. 2: SM contribution to the denominator in the differential forward-backward asymmetry as defined in 



( 18 1. CTEQ5M parton distribution functions were used. 



and is shown in Fig. (|2]). We note here that the falling SM contribution alone is not enough to give 
as steep a rise in the asymmetry as a function of s as is observed at CDF. The rise can steepen 
through a combination of the following factors: (1) | (T^^~ rises as a function of s and/or (2) 
i^ua'^ is comparable to the SM contribution and decreases as a function of s. However, if the 
majority of the steepness were to come from mechanism (2), the total cross-section especially at 
low invariant mass would have to be comparable to the SM cross-section; this is hard to do without 
running into constraints on the total differential cross-section. Thus a significant contribution must 



NP- 



come from a^^ 

There are seven kinds of terms that show up in a general cross-section involving t-channel 
mediators, including its interference with the SM: 



uf + smf if + smf uf if 



''M "-Af "-M 






(21) 



with t ^ n for or ^-channel diquarks. We examine these contributions term by term to determine 
which types can successfully generate a large contribution. In particular, there must be a large 
contribution to the asymmetry with a very modest contribution to the total cross-section. That is 
to say simply that the odd contribution must be large in comparison to the even contribution. 

We examine this in detail in Fig. (|3]) for different types of effective vertices. The salient points 
to take away from the figures are: (1) Scalars have odd contributions comparable to vectors only 
in the higher mediator mass range. (2) As a function of energy, the magnitude of the odd term 
for a given contribution is never greater than the magnitude of the even term, though some terms 
obtain much closer to equal magnitudes than others. Thus in order to best succeed in generating 
a sizable positive asymmetry while not destroying the invariant mass distribution, an ideal model 
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FIG. 3: Terms contributing to cross-sections with t or n channel mediators. Solid lines indicate the odd 
contribution and dotted the even contribution, integrated over cos Q. The top plots include terms from the 
interference term, and the bottom plots from the NP squared term. For diquarks, u -H- t, which flips the sign 
of the odd contribution and leaves the even contribution the same. The letters in square brackets indicate 
whether the term appears for scalar [S], vector [V], or both [S,V] mediators. 

will involve destructive interference betvv'een the even parts of the SM-NP interference and NP 
squared terms of the amplitude, and minimal or constructive interference between the odd part of 
the SM-NP interference and NP squared terms. ^ By inspection, none of the scalars (t or n-channel) 
can satisfy this condition. Scalar diquarks have some success in generating a substantial asymme- 
try in an intermediate mass range where the squared term contributes the dominant positive odd 
contribution. For the triplet, the interference term gives a negative odd and even contribution (so it 
helps to lower the cross-section but also lowers the numerator) while for the sextet the interference 
term enters with a minus sign and so gives a positive odd and even contribution (so it increases the 
numerator but also the cross- section). 

For vectors there are more terms in play, so the story is a bit more complicated. To show the 
effects on the total asymmetry, we plot the total asymmetry (and, when relevant, cross-sections) 
for all t- and w-channel mediator color representations and spin combinations in Figs. Q-Q. 



That some amount of destructive interference is favored by the data was noted in B6ll . 
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We show three benchmark mediator masses. Fig. (|4]) shows the scalar models that succeed in 
generating a positive asymmetry, though in general for perturbative couplings it is not a large 
positive asymmetry; Fig. ([5]) shows the same for the vector mediator case, and it is seen that the 
contribution to the total asymmetry can be large for all mass ranges. Lastly, we show in Fig. ([6]) 
the mediators that fail to produce a positive asymmetry larger than 5%. These include the scalar 
color octet and vector triplet and sextet. 

One might wonder whether the asymmetry induced by scalars could be enhanced by adding 
another scalar with s-channel couplings to uu and ti. This is predicted, for example, by the flavor 
triplet models. Interference between a t-channel scalar with mass rui and an s-channel scalar with 
mass mo would give rise to terms of the form ,, ^K.^: and ,. '^'^i.i . These contributions, assuming 
mi = 1712, are shown in Fig. (|7]). For mediators lighter than the top quark, the odd contribution 
has the same sign as even for both terms, and it is hard to see how these contributions can enhance 
the asymmetry while not increasing the total cross-section to unacceptable levels. For mediators 
heavier than the top quark, odd and even contributions for the s4/(sm^a4^) have the opposite sign — 
the odd contribution is positive for energies below the mediator mass and negative above. This 
interference could have interesting implications for models involving both t-channel and s-channel 
scalars of intermediate mass. Diquarks with s-channel interactions would not contribute to the ti 
cross-section or ApB- 

Lastly, we briefly discuss s-channel mediators, which can give rise to a large asymmetry for 
an appropriate choice of couplings and masses. If the asymmetry is generated from s-channel NP 
interactions, then the cos 9 dependence of the NP cross-section is a simple quadratic polynomial. 
The axigluon originally proposed in Ref. [fTSl [T9l supposed a heavy axigluon to evade dijet and 
ti resonance searches that strongly constrain the state with masses below ~ 2 TeV. However, 
recently different regions of the axigluon parameter space have been explored. For example, a 750 
GeV state was considered in [35l, with the dijet constraints evaded by making the coupling to the 
top quark much larger than to the up quarks. A 400 GeV state was considered in ||3T| , and the 
ti resonance search constraints evaded by making the state sufficiently broad. Lastly, if a vector 
with diagonal axial couplings to top and up has a mass slightly lighter than the top mass, then it 
will not show up as a resonance in the ti spectrum. We refer the reader to these references for 
details, though we include the axigluons in our scans of parameter space in the next section for 
completeness. 
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FIG. 4: Spin-0 mediators. Left-hand plots show the differential asymmetry for various couplings given 
a 150 GeV, 400 GeV or 800 GeV mediator A line is drawn at Vs = 450 GeV to highlight the value 
of the asymmetry at the lower end of the the CDF analysis higher invariant mass bin. Due to the rapidly 
falling PDFs, the high invariant mass bin asymmetry will be given roughly by the value of the differential 
asymmetry at 450 GeV. The right-hand plots show contributions to the parton level uu — )■ tt cross-section as 
a function of \/I (dotted lines), and to the odd parton level cross-section (forward - backward), normalized 
by 327rs//3 to make a dimensionless quantity. The effective Standard Model contribution as defined in 
( |20] ) is shown as a black dotted line. Contributions to the total differential cross-section, dai{s)/ds, can be 
obtained from the dotted contributions by multiplying by the factor §2^- ^^^^ ^1^- < P^ ^i^d ( [T9] )) 
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FIG. 5: Spin-1 mediators. Left-hand plots show the differential asymmetry for various couplings given 
a 150 GeV, 400 GeV or 800 GeV mediator. A line is drawn at \/I = 450 GeV to highlight the value 
of the asymmetry at the lower end of the the CDF analysis higher invariant mass bin. Due to the rapidly 
falling PDFs, the high invariant mass bin asymmetry will be given roughly by the value of the differential 
asymmetry at 450 GeV. Right-hand plots show contributions to the parton level uu — )• tt cross-section as a 



function of v s, as in Fig. |4j The Standard Model contribution as defined in ( |20| ) is shown as a black dotted 
line. 
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FIG. 6: The asymmetry for representations that cannot produce a positive asymmetry of more than a few 
percent. 
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FIG. 7: Terms contributing to s-channel scalar / i-channel scalar interference cross-section. Solid lines are 
the odd contributions and dashed are the even contributions, integrated over cos Q. Here we assume a narrow 
width. 

ni. COMPREHENSIVE SEARCH FOR MODELS AT THE TEVATRON 

To augment the conclusions of the previous section, we carry out a comprehensive representa- 
tive scan of models using MadGraph[l70l|; the details of our procedure are discussed in Appendix 
|B} We scan over s, t, w-channel models, characterized by a single new mediator of given spin and 
color representation (|2]); we scan over all such models that can produce a positive asymmetry of 
more than a few percent while remaining (somewhat) perturbative (coupling < 6) and contributing 
less than order 50% to the total cross-section in the mass range 200 GeV - 2 TeV."' The models 



scanned are summarized in Table III We choose representative models that generate the largest 



asymmetry. For t-channel models we focus on mediators connecting up to top, both because they 
generate a large asymmetry, and also because a light neutral state runs into few constraints. The 
color singlet and triplet are our representative scalar models, though neither is successful in gen- 
erating a large asymmetry, as we detailed earlier. Also note that the singlet scalar is part of an 
electroweak doublet, though we choose to couple this scalar to t^ — ur so that only one state is 
operative for the forward-backward asymmetry. The charged component of the SU{2) mediator 
multiplet will contribute to bb plus jet events at the LHC, but this will be easily overwhelmed 
by the background. For the t-channel flavor-violating models, we consider both a color singlet 
vector (CIV) and octet vector (CSV) that couples only to right-handed states. We also consider a 
flavor octet, color singlet vector (F8C1V) that couples to Un^y^Un, where now the up quarks are 

'° We neglect models with mass below the top {e.g. Il50l ): in general, these models will tend to rather severely 
overproduce the total cross-section and number of additional jets at the LHC and/or lead to large contributions to 
single top production, depending on the details of the mediator decay channels. 



15 



Model 



Spin Color SU{2)y Flavor s-, t-, u-7 Comments and References 



CIS 







1 



'1/2 



1 



C3S 





3 


I4/3 


1 


u 


CIV 


1 


1 


lo 


1 


t 


CSV 


1 


8 


lo 


1 


t 


F8C1V 


1 


1 


lo 


8 


t, s 



Only very moderate asymmetries achiev- 
able 0{'> 10%). Low mass (mj\/ ~ mt) 
states do slightly better ISTI . 
a.k.a. triplet diquark. q = 4/3 
a.k.a. Z' or W 1501. 



schanC8V(A,R) 1 



schanC8V T 



Flavor breaking only through up Yukawa 

m. 

a.k.a. axigluon or coloron. For 27nt < 
ruM ^ 2TeV, very broad width required 
to avoid tt resonance searches lfT9l[35ll56l . 
~ 400 GeV broad resonance via additional 
scalars. Universal quark couplings [31]. 



TABLE IIL Summary of models scanned. All t- or u- channel states are taken to be non-self-conjugate. 

in an octet of SU{3)uj^. Lastly, the s-channel axiglue type models are considered, both in flavor 
universal [1311 and non-universal [[T9ll35ll varieties. 

The results of this scan for the Tevatron are shown in Figs. ([8]) - (fTO]). The coupling conventions 



in the figures are as follows. The t-channel scalars, as well as CI V and CSV, models are labeled by 
their coupling to RH quarks, with gi = 0- The flavor symmetric F8C1V model has an additional 
parameter r] that controls the flavor breaking coupling to the top quarks such that couplings to 
top-quarks have couplings ^ gn + 2r]m'^/v^, with v = 246 GeV. The coupling conventions for 
s-channel models are more complicated. The couplings in schanCSVF and schanCSVA are purely 
axial {gR = —gi), with the former only being flavor universal. The schanCSVF model has an 
independent width parameter [|31ll . which was scanned over to find models with maximally large 
asymmetries per unit production cross-section. schanCSVR has non-universal couplings to right- 
handed quarks [|35l . 

We apply cuts on the simulated sample and fully reconstruct tops as described in Appendix |B] 
to mimmic the analysis in fT\. More specifically, for Tevatron events we apply the following sets 
of cuts: 
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signal level Am at Tevatron, cr = £riev''Accep!ai] 




FIG. 8 : Scatter plots depicting simulated signal level {ApB {rrif^ < 450 GeV) , Aps [m^^ > 450 GeV) , a^i x 



acceptance} at Tevatron CM energy for t-channel flavor-changing scalar models listed in Table III The 
models are labeled by the mass of the mediator, the coupling to right-handed quarks, and the total Tevatron 
production cross-section times acceptance. The cross-sections are compared against the SM cross-section 
times acceptance which yields 0.252 pb at the LO; the color scales for the models indicate the deviation 
from the SM cross-section, as indicated by the legend at the bottom. The curves indicate constant x^ for 
a given cross-section, as defined in Eq. ^27) . Contours for four cross-section values (cyan, blue, green, 
purple) are shown for x^/d.o.f. = 1 and 2. A single (cyan) x^/d.o.f. = 3 contour with SM cross-section is 
shown. Model points of a given color should be compared to x^ contours of the same color. 
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signal level Am at Tevatron, cr = £riev''Accep!ai] 




FIG. 9: Scatter plots depicting simulated signal level {^^^^(Tnjf < 450 GeV) , ^i^^ ("i^tf > 450 GeV) , cr^j x 



Acceptance} at Tevatron CM energy for t-channel flavor-changing vector models listed in Table III The 
models are labeled by the mass of the mediator, the coupling, and the total Tevatron production cross- 
section times acceptance. The coupling conventions are discussed in detail in the text. The cross-sections 
are compared against the SM cross-section times acceptance which yields 0.252 pb at the LO; the color 
scales for the models indicate the deviation from the SM cross-section, as indicated by the legend at the 
bottom. The curves indicate constant x^ for ^ given cross-section, as defined in Eq. (|22]). Contours for 
four cross-section values (cyan, blue, green, purple) are shown for x^/d.o.f. = 1 and 2. A single (cyan) 
X^/d.o.f. = 3 contour with SM cross-section is shown. Model points of a given color should be compared 
to x^ contours of the same color. 
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signal level Am at Tevatron, cr = £riev''Accep!ai] 




FIG. 10: Scatter plots depicting simulated signal level {ApsiiTT'tt < 450 G&N),AFB{mtt > 



450 GeV) , ati x Acceptance} at Tevatron CM energy for axigluon models listed in Table III The mod- 
els are labeled by the mass of the mediator, the coupling, and the total Tevatron production cross-section 
times acceptance. The coupling conventions are discussed in detail in the text. The cross-sections are com- 
pared against the SM cross-section times acceptance which yields 0.252 pb at the LO; the color scales for 
the models indicate the deviation from the SM cross-section, as indicated by the legend at the bottom. The 
curves indicate constant x^ for a given cross-section, as defined in Eq. (|22]). Contours for four cross-section 
values (cyan, blue, green, purple) are shown for x^/d.o.f. = 1 and 2. A single (cyan) x^/d.o.f. = 3 contour 
with SM cross-section is shown. Model points of a given color should be compared to x^ contours of the 
same color. 
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• Exactly one electron or muon with pr > 20 GeV and \ri\ < 1.0. 

• At least four jets with pr > 20 GeV and \r]\ < 2.0, with at least one of the jets having a 
b-tag. 

• Ef'' > 20 GeV. 

We reconstruct tops as described in [61 J, doing a likelihood analysis on the lepton and jet kine- 
matics to the tt hypothesis, using the algorithm described in our previous paper [l63l|. The cone 
jet algorithm was used for Tevatron events. Jet energy scale corrections were carried out via a 
procedure described in Appendix |B] 

We choose to show results after detector simulation (at the signal level) because, as discussed 



in [16111 . unfolding of data to the parton level is model dependent. In Figs. ([8pO|) the axes give the 
signal level ApB with rritt < 450 GeV and mu > 450 GeV. The ellipses encircle the best fit points 
to the CDF signal level semileptonic ti Aps with concentric ellipse giving x^/d-o.f. = 1, 2, 3, 
with the constraints from the total cross-section times acceptance being taken into account via 

2 ^ I ^FB ~ ^FB obs j , I ^FB ~ ^FB obs j , / ^tt - C^ff ,SM \ .^2) 



We use 






cr^™,,<450 = VO.O392 + 0.0172 + 0.0252 + 0.0152 (23) 

^FB 

a^™^,>450 = VO.O532 + 0.0322 + 0.0252 + 0.0432 (24) 

^FB 

^"JL. = 0. 10 (LHC7) or 0. 15 (TEV) (25) 



0'tt,SM 

for the error estimates and 

^Fr''°obs = -0-022 (26) 

^7f\b. = 0-266 (27) 

(rti,sM = 8.18 pb (LHC7) or 0.252 pb (TEV) (28) 

as the central values. Note that the last value is the central LO SM cross-section times acceptance, 
given the cuts for Tevatron and LHC7 outlined in this section and the next. The central values 
for ApB are the background- subtracted signal level values from [7J and the SM values for the 
cross-section times acceptance are taken from our simulations of 5 million events. The first two 
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contributions to the ApB errors are from experiment, the third for the typical statistical error from 
our finite-sized simulated data samples, and the last is to account for possible NLO corrections: 
we take this contribution to be of the same size as the NLO SM asymmetry. For the cross-section 
error, we take 15% errors for the Tevatron and 10% for LHC. 10% is roughly the current ex- 
perimental error for the Tevatron measurements, and we add a ~ 5% uncertainty due to the top 
mass and theory uncertainties in the NLO corrections. For LHC, the statistics on the cross-section 
measurement should lead to smaller error bars and we take this into account with a smaller LHC 
error of 10%. A value x^ ~ 3 indicates a good fit to data. For the Standard Model with ApB 
given by the NLO prediction and cross-section by our LO simulations, x^/3 = 2.8. Since we take 
the central value for the cross-section to be the SM LO value, this value is somewhat artificially 
low. These error estimates should be taken as rules of thumb to guide the eye in our figures for 
comparing SM against NP, rather than as hard and fast quantitative error budgets. 

We discuss the scalar models first. As can be seen from Fig. ([8]), the triplet scalars generally 
produce larger asymmetries than singlet scalars, which generally cannot produce a larger asymme- 
try than the SM. This can be qualified if the singlet scalars are lighter than the top mass, in which 
case signal level asymmetries as large as 10% for mu > 450 GeV can be achieved (though this is 
well below what is observed). This in agreement with the parton level results of [51], taking into 
account the factor ~ 2 washout translating from parton level to signal level. The triplet scalars 
seem to reproduce the total asymmetry and cross-section very well. However, it was shown in I16f]| 
that these models seriously overproduce the invariant mass distribution at large invariant mass. We 
refer the reader to [|6T1l for details. 

Next we discuss the t-channel vector mediators in Fig. (|9]). As expected from the results in [|6TI . 
the color singlet vector is most successful in reproducing the asymmetry at high invariant mass and 
satisfying the cross-section constraints. Due to details in the form of the matrix element, the color 
octet is less successful. The flavor universal octet can produce large asymmetries, but these also 
tend to come with fairly large contributions to the total cross-section, due to the presence of both 
s and t-channel mediators. 

Lastly, we discuss the s-channel states in Fig. ( fTO] ). The wide, low mass axiglue models, 
schanCSVr, in general are most successful at producing a large asymmetry with small contri- 
bution to the total cross-section. The light axigluon models with couplings to right-handed quarks 
and masses in the 700-900 GeV range (schanCSVR) [13511 do not produce a large asymmetry on 
the other hand; in most cases it is not larger than the SM asymmetry. Heavy axigluon models can 
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succeed with a large enough coupling to light quarks, but these risk being ruled out shortly by 
LHC dijet and tt resonance searches. 

With these results in hand, we now turn to examining the implications of models that are capa- 
ble of satisfying the Tevatron constraints on top analyses at LHC7. For each class of models, and 
a selection of mediator masses between 200 GeV and 2 TeV, we take models with the lowest x^ 



as defined by the statistic in Eq. (22 1. 5 million events are generated for each of these benchmark 
models to gain enough statistics at the high invariant mass, via the procedure in AppedixJB] Our 
benchmark models are not an exhaustive set of model choices, but they are indicative of the types 



of models that can generate top Aps ■ The choice of models is shown in Table IV It gives the mass 
and coupling of the model, the LO cross-section at Tevatron and LHC along with the acceptance 
A., the signal and parton level ApB in the low and high invariant mass bins, along with the total 
asymmetry, and the x^ ^t Tevatron and LHC, using the statistic discussed in the text. These results 
are to be compared against the SM, shown in Table [V] 

Note that the models with the lowest x^ tend to universally underproduce the total asymmetry. 
The reason is that the models with the largest ApB also tend to overproduce the total cross-section 
rather seriously, so that the x^ prefers to take a hit on the asymmetry (which has la errors of 
~ 8% at the signal level) in lieu of a large ti production cross-section. The models that are the 
least successful at producing a large asymmetry with minimal impact on the total cross-section 
are: CSV, CIS, schanCSVR. These models are generally able to produce little more than the SM 
asymmetry for ApB with rritt > 450 GeV, and should not be considered as viable models for Aps- 

Before moving on to the LHC analysis, we check the Tevatron invariant mass distributions 
for the classes of models that we examine more carefully. As we learned in |6T1, acceptance 
effects can be important in bringing NP models into agreement with the Tevatron invariant mass 



distributions. We show the Tevatron invariant mass distributions in Figs. (|TT])-([T2|), for comparison 
to the LHC results we discuss next. 

IV. IMPLICATIONS FOR TOP PHYSICS AT THE LHC 

For the LHC benchmark points analysis, we generate 5 million events for each model, as we 
did for the Tevatron analysis. We also modified the PGS code to implement the anti-kp algorithm 
[TTTII to mimmic ATLAS as detailed in Appendix |Bj In the following, closest attention should be 
paid to the CIV, F8C1V, CSS, schanCSVF and schanCSVA models, as these, among the models 
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signal level parton level 

parameters axev (pb), A. ^lhc (pb), A. A<f^, A>g°, Af^' A<f ", A>f °, Af^' xL^ xjuc 
m, QR CIV 



200.. 


,0.7 


6.3, 


0.037 


146, 0.068 


-0.03,0.15,0.06 


0.01,0.39,0.2 


0.8, 


2.3 


400.. 


,1.3 


7.1, 


0.038 


154, 0.073 


0.01,0.25,0.15 


0.08, 0.55, 0.35 


0.2, 


4.7 


600.. 


,1.5 


5.3, 


0.039 


126, 0.072 


-0.04,0.15,0.06 


-0.03,0.25,0.1 


1.2, 


1.2 


800.. 


,2.1 


5.8, 


0.039 


129, 0.073 


-0.03,0.18,0.09 


-0.01,0.36,0.18 


0.5, 


1.1 



m, QR C8V 



400., 0.75 6.8,0.041 130,0.072 0.01,0.08,0.04 0.03,0.1,0.06 2.1,2.7 

800., 1.4 6.8,0.04 120,0.072 -0.01,0.08,0.03 0.03,0.1,0.06 1.9,2. 



m,g,ri F8C1V 



200., 0.5, 1. 6.5,0.037 148,0.067 0.05,0.14,0.09 0.03,0.4,0.21 1.4,2.8 

400., 0.5,0. 9.4,0.04 125,0.069 0.08, 0., 0.05 0.23,-0.02,0.17 8.5,5.1 

600., 0.5, 3. 6., 0.041 128,0.071 -0.03,0.15,0.07 -0.05,0.31,0.14 0.7,1.1 



800., 0.5, 1. 


6., 0.041 


115,0.072 


-0.03,0.01,-0.01 


0., 0.03, 0.01 


3.5, 3.5 


m,gR 






CIS 






200., 1.5 


5.7, 0.042 


119,0.072 


0.01,0.04,0.03 


0., 0.06, 0.02 


2.9, 3. 


m,gR 






C3S 







400., 2.95 8.6,0.033 165,0.074 0., 0.17, 0.11 0.2,0.22,0.21 0.8,8.4 

600., 3.4 6.7,0.043 133,0.075 0., 0.14, 0.08 0.05,0.23,0.14 1.2,2.6 

800., 4.15 6.6,0.042 128,0.075 -0.01,0.15,0.08 0.03,0.27,0.15 0.9,1.8 



m, gn, r/m(%) 






schanC8Vr 






420., 0.45, 18 
440., 0.45, 13 


6.7, 0.04 
6.9, 0.039 


116,0.072 
118,0.07 


-0.03,0.15,0.05 
-0.03,0.12,0.04 


-0.03,0.3,0.1 
-0.11,0.34,0.06 


0.8,0.8 
1.1,1.1 


^' 9% 9r 






schanC8VA 







2000., -1., 5. 6.4,0.04 117,0.072 0.01,0.16,0.08 

2400., -3.6, 3.6 6.5,0.039 119,0.072 0., 0.14, 0.07 

m, gj^, gj^ schanC8VR 

700., -0.05, 4.5 6.7,0.04 116,0.07 0., 0.06, 0.02 

850., -0.08, 6. 6.7,0.039 117,0.072 0.04,0.08,0.06 

TABLE IV: A representative set of models chosen for LHC analysis. Acceptance is labeled "A.". 
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0.06,0.17,0.1 


0.7, 


0.8 


0.07,0.21,0.13 


1., 


1. 


0.02, 0.07, 0.04 


2.4, 


2.4 


0.02, 0.08,0.04 


2.2, 


2.2 



LO SM cross-section, NLO ApB 



signal level parton level 

cTTev (pb), A. cTLHC (pb), A- A<f^, A>f^, A^^ ^<g°, A>g°, Af^^ xjy^ Xlnc 
6.3,0.04 115,0.071 0.015,0.043,0.024 0.040,0.088,0.058 2.8,2.8 

TABLE V: The SM LO cross-section at Tevatron and LHC along with the acceptance. A., the signal and 
parton level ApB in the low and high invariant mass bins, along with the total asymmetry, and the x^ at 
Tevatron and LHC, using the statistic discussed in the text. 
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FIG. 11: Tevatron invariant mass distributions, on both linear and log scales, for our benchmark models 
choices. The SM is shown in the yellow band, with statistical errors for 5.3 fb^^ of data. 
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FIG. 12: Tevatron invariant mass distributions, on both linear and log scales, for our benchmark models 
choices. The SM is shown in the yellow band, with statistical errors for 5.3 fb~^ of data. 

in the literature we have considered, are able to generate the top ApB to a reasonable degree. 
The variables that we focus on at the LHC are: 

• Total cross-section. The chief uncertainties here come from NLO corrections from both the 
SM and NP, and the uncertainty in the top mass; 

• Invariant mass distribution. Here again NLO corrections will play an important role; 

• Number of additional jets. In t-channel models, single production of the mediator in con- 
junction with the top is an important process at the LHC. A gluon and light quark in the 
initial state will exchange a top in the t-channel and produce a top along with a mediator as 
in Fig. ([!]). The mediator will prefer to decay to a top and another jet, leading to a potential 
enrichment of events with an extra jet. The direct search for the top-jet resonance as a signa- 
ture for these models was studied in [|63l . but its presence may be known through counting 
the number of additional jets in tt events. 

• Rapidity distribution of the lepton. Especially for models with a t-channel resonance, the 
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leptons may be produced in a more forward direction at the high invariant mass. On the 
other hand, single mediator production leading to ti + jets events can lead to more central 
leptons. 

We follow the cuts discussed in the 200 pb~^ ATLAS semileptonic top analysis flSl. We require: 

• exactly one electron with pt > 25 GeV and \r]\ < 2.5, or exactly one muon with px > 
20 GeV and \t]\ < 2.5; 

• at least four jets with pt > 25 GeV and \r]\ < 2.5, one of which must be 6-tagged; 

• if the lepton is an electron, we require EJp^'^'^ > 35 GeV and the transverse mass of the lepton 
and E^^'^^ be greater than 25 GeV; if the lepton is a muon, we require i?^*** > 20 GeV and 
the transverse mass of the lepton with Elp^^^, plus the E^^'^^, be greater than 60 GeV; 

• jets within AR < 0.2 are removed so as to avoid double-counting of electrons as jets. 

In addition, ATLAS demands isolation cuts; since we do clustering in PGS before placing the 
cuts, we do not apply them, m^f is re-constructed in the same way as ATLAS, carried out without 
a full top reconstruction. The neutrino momentum is found assuming the W mass and massless 
neutrino conditions. For some events there is no positive energy solution, in which case the event 
is discarded. According to the ATLAS analysis, we take the longitudinal neutrino momentum to 
be the real part of the mass constraint solution in the case of imaginary solutions and we take the 
solution with smallest absolute value if there are two solutions. 

The first and simplest measure is the top forward-backward asymmetry versus the total produc- 
tion cross-section at the LHC. There is a trade-off between models with a large enough coupling to 
produce the observed forward-backward asymmetry, while simultaneously having a small enough 
coupling that single mediator production at the LHC does not lead to a large contribution to the 
ti cross-section. However, given that the higher mass models in particular have large couplings, 
one expects the next-to-leading-order (NLO) corrections to play a significant role in both the total 
cross-section and invariant mass distributions. Given the gross-overproduction of some models 
of the total cross-section, some may, however, be reasonably eliminated. This can be seen in 



Figs. (13 15), where we plot the Tevatron ApB in low and high invariant mass bins, with total 



production cross-section at LHC times acceptance indicated by color. We again compare LO 



MadGraph results against the LO SM cross-section times acceptance ( |28] ). We find a total LO 
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matched SM tt + or 1 jets cross-section of 1 15 pb for nit = 172 GeV. Note that there is a large 
K-factor of ~ 1.6 expected at LHC7 which enters to match the total cross-section observed (of 
about 180 pb) against the LO contribution. The LO cross-section times acceptance is 8.178 pb. 



A couple of features in Figs. 13 15 1 in particular are worth highlighting, t-channel models at 
the LHC overproduce the total cross-section much more than at the Tevatron. This is because 
single production of the t-channel mediators gives rise to a significant contribution to the total 
cross-section. This effect is more important for lighter mediators, so that light mediators become 
much more disfavored at the LHC. 

Note that this brings in a significant tension for t-channel models of Aps between the con- 
straints from Tevatron and the LHC. The Tevatron ti invariant mass distribution tended to favor 
lighter mass mediators because they lead to less distortion of the invariant mass distribution at high 
invariant mass ||6TI . while LHC favors heavier mediators because they lead to less distortion in the 
total cross-section. 

We next consider the effect on the invariant mass distribution for our benchmark models. The 



results are shown in Figs. ( [16] ), ( [T7| ). The effect of the NP on the shape of the invariant mass 
distribution is very different at the LHC than at the Tevatron. At the Tevatron, the effects of the 
new mediator become most pronounced at the high invariant mass — for the t-channel models in 
particular. At the LHC, this effect is not present, because most of the impact of the new mediators 
on the cross-section is simply single mediator production. The s-channel models with a sufficiently 
broad width have little impact on the invariant mass distribution. Thus it appears that most of the 
constraint on the new models comes simply from the total cross-section measurement. We also 
note that acceptance effects explored in []6T1| are not as important at LHC as at Tevatron, both 
because NP ti events at LHC are more central, and because the rapidity coverage for leptons at 
ATLAS and CMS is better than at Tevatron. 

Perhaps the leading discriminant is simply the number of additional jets in the event, shown 



in Figs. ([18]), ( [T9| ). While the overall production cross- section may be somewhat uncertain due to 
NLO corrections, leading to an uncertainty in the overall normalization of the NP curves, there 
is a significant difference in the ratio of the number of events with one extra jet to the number 
of events with no extra jets. In fact, all of the t-channel models that generate a large asymmetry 
significantly overproduce the number of events with one additional jet. One might wonder whether 
this effect could be reduced by allowing a significant branching fraction to light quarks; however 
in this case these events will contribute significantly to the single top analyses, which already 
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signal level ApB at Tevatron, cr = trLHC''Acceptar 




FIG. 13: Scatter plot depicting simulated signaIIeveI{Ai7'B(m4f < 450 GeV) , Ai^s (Tn^f > 450GeV), u^fX 



acceptance} for t-channel flavor-changing scalar models listed in Table III The models are labeled by the 
mass of the mediator, the coupling to right-handed quarks, and the total LHC production cross-section times 
acceptance. The cross-sections are compared against the SM cross-section times acceptance which yields 
8.178 pb at the LO; the color scales for the models indicate the deviation from the SM cross-section, as 
indicated by the legend at the bottom. The curves indicate constant x^ for a given cross-section, as defined 
in Eq. ([22]). Contours for four cross-section values (cyan, blue, green, purple) are shown for x^/d.o.f. = 
I and 2. A single (cyan) x^/d.o.f. = 3 contour with SM cross-section is shown. Model points of a given 
color should be compared to x^ contours of the same color. 
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signal level ApB at Tevatron, cr = trLHC''Acceptar 




FIG. 14: Scatter plot depicting simulated signaIIeveI{Ai7'B(m4f < 450 GeV) , Air^ ("^if > 450GeV), cr^fX 



Acceptance} for t-channel flavor-changing vector models models listed in Table III The models are labeled 
by the mass of the mediator, the coupling, and the total LHC production cross-section times acceptance. 
The coupling conventions are discussed in detail in the text. The cross-sections are compared against 
the SM cross-section times acceptance which yields 8.178 pb at the LO; the color scales for the models 
indicate the deviation from the SM cross-section, as indicated by the legend at the bottom. A single (cyan) 
X^/d.o.f. = 3 contour with SM cross-section is shown. The curves indicate constant x^ for a given cross- 
section, as defined in Eq. ( [22| ). Contours for four cross-section values (cyan, blue, green, purple) are shown 
for x^/d.o.f . = 1 and 2. Model points of a given color should be compared to x^ contours of the same color. 
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signal level ApB at Tevatron, cr = trLHC''Acceptar 




FIG. 15: Scatter plot depicting simulated signaIIeveI{Ai7'B(m4f < 450 GeV) , Ai^^ ("i^if > 450GeV), u^fX 



Acceptance} for axigluon models listed in Table III The models are labeled by the mass of the mediator, 
the coupling, and the total LHC production cross-section times acceptance. The coupling conventions 
are discussed in detail in the text. The cross-sections are compared against the SM cross-section times 
acceptance which yields 8.178 pb at the LO; the color scales for the models indicate the deviation from the 
SM cross-section, as indicated by the legend at the bottom. A single (cyan) x^/d.o.f. = 3 contour with SM 
cross-section is shown. The curves indicate constant x^ for ^ given cross-section, as defined in Eq. (22 1. 
Contours for four cross-section values (cyan, blue, green, purple) are shown for x^/d.o.f. = 1 and 2. Model 
points of a given color should be compared to x^ contours of the same color. 
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FIG. 16: CIS, C3S, CIV, F8C1V models differential cross-section times acceptance (on both linear and 
log scales) at LHC7 versus reconstructed m^f, as compared to SM LO expectation, with ztlo" yellow bands 
corresponding to statistical error given Ifb^^. Models shown are those with the lowest x^ for a given mass 
as defined in Eq. ([22]), except for the 600 and 800 GeV CIV models, which were chosen to have the lowest 
X^ and be within 10% of the SM cross-section times acceptance at Tevatron. 

with only 200 pb~^ of data have an uncertainty of only 40 pb [TtTII . In the case of a significant 
branching fraction of the mediator to light quarks, single mediator production will easily contribute 
a significant fraction of this cross-section, with even more severe constraints arising in the high Ht 
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FIG. 17: CSV and schanCS models differential cross-section times acceptance (on both linear and log 
scales) at LHC7 versus reconstructed m^j, as compared to SM LO expectation, with itlcr yellow bands 
corresponding to statistical error given Ifb^^. Models shown are those with the lowest x^ as defined in 
Eq. ([22]l. 



tail of the distribution, as pointed out in [48] . For reference in these figures we have also shown 
the rapidity distribution of the leptons; single mediator production results in a more central lepton 
rapidity distribution. 



V. CONCLUSIONS 

We have carried out a comprehensive analysis of NP models for top ApB utilizing both Tevatron 
and the prospective LHC7 constraints with 1 fb^^. We considered effective vertices for all pos- 
sible spin, color and flavor representations connecting top quarks with up quarks. We were able 
to show on general grounds why scalar mediated models have difficulty reproducing the observed 
asymmetry. We revisited the Tevatron signal level invariant mass distributions, as investigated in 
our earlier paper [61 1. We found that the prospective LHC constraints on the total cross-section of- 
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FIG. 1 8 : Number of additional jets and lepton differential rapidity distribution of C 1 V, F8C IV and schanCS 
models at LHC7, as compared to SM LO expectation, with zbla yellow bands corresponding to statistical 
error given Ifb^^. Models shown are those with the lowest x^ ^s defined in Eq. ( [22] ). 

far complimentary constraints to the Tevatron invariant mass distribution. In the case of t-channel 
mediators, the LHC total cross-section places a strong constraint on light mediators, while the 
Tevatron invariant mass distributions place strong constraints on heavy mediators that are able to 
produce the asymmetry. The vanilla t-channel models thus seem disfavored at present. Heavy, 
narrow axigluons (with masses ~ 2 TeV) are currently becoming more tightly constrained with 
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FIG. 19: Number of additional jets and lepton differential rapidity distribution of schanCS, CSV, CIS 
and C3S models at LHC7, as compared to SM LO expectation, with itlo" yellow bands corresponding to 
statistical error given Ifb^^. Models shown are those with the lowest x^ ^s defined in Eq. ([22]). 



the recent LHC7 top results. A 400 GeV axigluon with large width and universal couplings to 
quarks appears at present to evade all existing constraints. 

The LHC is rapidly closing the window on viable models for the top forward-backward asym- 
metry. More non-generic features, such as large widths as in the light axigluon discussed here, 
will be necessary to make viable models consistent with both the Tevatron top Aps and LHC top 
observables. 
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Note added 1: After the first version of this manuscript was completed, a new result on top 
ApB from the DO collaboration appeared 11721 . in which ApB in the high tt invariant mass bin is 
significantly lower than that of the CDF result. As a result, the concentric x^ contour ellipses in 
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Figs. ([8 10), ( 13 15 1, and (20 22) will move down and to the right when the CDF and DO results 



are combined, so that many model points in danger with CDF alone will have a significantly lower 
X^. As a result, the best model point may change. We leave the analysis to a future publication. 
Also after the first version of this manuscript was completed, new results on tt cross-sections at 



LHC7 were released at the 201 1 International Europhysics Conference on High Energy Physics 



No significant deviation from the SM expectation was measured 11731 . 

Note added 2: Because of the rapidly changing experimental results, we will periodically up- 
date our results on a web page. This website will also provide some figures not included in this 
paper. See http://susy.physics.lsa.umich.edu/TopPhysics, 

Appendix A: Parton Level Asymmetries 



As a complement to the Tevatron signal level asymmetries shown in Figs. rt8])-( 10 1, we show the 



parton level asymmetries, so that theorists can easily map signal level onto parton level for a broad 



range of models. These are shown in Figs. (|20|)-<[22]) for the same model points as in Figs. (|8|)-(|T0]). 
Note in comparing the signal and parton level plots that a number of points are deleted in the parton 
level plot in cases where they cluster strongly around the SM point and become indistinguishable. 



For these parton level plots the x statistic used to draw contours is defined in Eq. (22 1, but with 



a^,j<45o = VO.1462 + 0.0472 + 0.0052 + 0.0402 (Al) 



■^FB 



a^^,i>Aw = VO.IOI2 + 0.0492 + 0.0052 + 0.0882 (A2) 

^^^ = 0.15 (A3) 

0"tt,5Af 

for the error estimates, and with 

^^r''°obs = -0-116 (A4) 

AFf''\,s = 0-475 (A5) 

a„-sM = 6.27 (A6) 

for the central values. The central values for ApB are the parton level values from [7J and the 
SM values for the cross-section are taken from our simulations of 5 million events. The first two 
contributions to the ApB errors are from experiment, the third for the typical statistical error from 
our finite-sized simulated data samples, and the last is to account for possible NLO corrections: 
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we take this contribution to be of the same size as the NLO SM asymmetry. These error estimates 
should be taken as rules of thumb to guide the eye in our figures for comparing SM against NP, 
rather than as hard and fast quantitative error budgets. 

Appendix B: Event Generation 

In this appendix, we describe our event generation setup and strategies for data analysis 
presented in the main text of this paper. We employ FeynRules vl . 4 . 10 for model file 
generation [1741 . MadGraphS 1.3.3 for event generation [|70l . PYTHIA v6 for parton shower- 
ing and hadronization lTTSl . and a modified PGS4 for fast detector simulation [|76l . 

This work involves a large survey of different models and model parameters, and model- 
dependent acceptance in detection is an important issue in interpreting experimental observations. 
Thus fast detector simulation on a large number of events is necessary. Although there are criti- 
cisms on the credibility of fast detector simulation tools, fast detector simulation tools like PCS 4 
are indispensible for this paper. 

To obtain more realistic and reliable results, we tune the detector simulation and our analysis 
to the current experiments in such a way that performance is not harmed. For comparison of our 
results to data, we show NP models compared to the SM with the same analysis setup. Then we 
can draw conclusions on the status of NP models, since all experimental analyses are accompanied 
with their own SM simulation. In the following section, we summarize our considerations. 

1. Fast Detector Simulation and Object Reconstruction 

We simulate our model points given the specifications of the CDF detector at the Tevatron and 
from the ATLAS detector at the LHC. We use the default detector parameters for CDF and ATLAS 
given in the official distribution of PGS4. Some important detector parameters used in PCS 4 are 
summarized in Table IVTl 

While we have not modified the detector parameters, the PGS algorithm used for object re- 
construction is rather outdated and therefore can give rise to significantly different results. We 
summarize our changes in the following. 

a. Jet Reconstruction and Jet Energy Scale Correction: The official version of P G S supports 
two jet algorithms: the cone algorithm and A;T-jet algorithm. While the CDF analysis on ApB 
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FIG. 20: Scatter plots depicting simulated parton level {AFBitnti < 450 GeY), ApBifntt > 



450 GeV) , (Tjf } at Tevatron CM energy for t-channel flavor-changing scalar models listed in Table III The 



models are labeled by the mass of the mediator, the coupling to right-handed quarks, and the total Tevatron 
production cross-section. The cross-sections are compared against the SM cross-section which yields 6.3 
pb at the LO; the color scales for the models indicate the deviation from the SM cross-section, as indicated 
by the legend at the bottom. The curves indicate constant x^ for a given cross-section, as defined in Eq. 
( |22l ). Contours for four cross-section values (cyan, blue, green, purple) are shown for x^/d.o.f. = 1 and 
2. A single (cyan) x^/d-o.f. = 3 contour with SM cross-section is shown. Model points of a given color 
should be compared to x^ contours of the same color. 
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parton level .4|.|j at Tevatron. cr = crj^^. 




FIG. 21: Scatter plots depicting simulated parton level {ApB{mii < 450 GeY), ApBifnu > 



450 GeV) , a^t} at Tevatron CM energy for t-channel flavor-changing vector models listed in Table III 



The models are labeled by the mass of the mediator, the coupling, and the total Tevatron production cross- 
section. The coupling conventions are discussed in detail in the text. The cross-sections are compared 
against the SM cross-section which yields 6.3 pb at the LO; the color scales for the models indicate the de- 
viation from the SM cross-section, as indicated by the legend at the bottom. The curves indicate constant x^ 
for a given cross-section, as defined in Eq. (|22]l. Contours for four cross-section values (cyan, blue, green, 
purple) are shown for x^/d.o.f. = 1 and 2. A single (cyan) x^/d.o.f. = 3 contour with SM cross-section is 
shown. Model points of a given color should be compared to x^ contours of the same color. 
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parton level Ai.-u a! Tevatron. cr = crj^.^, 




FIG. 22: Scatter plots depicting simulated parton level {AFBitnti < 450 GeY), AFBifntt > 



450 GeV),(j4f} at Tevatron CM energy for axigluon models listed in Table III The models are labeled 



by the mass of the mediator, the coupling, and the total Tevatron production cross-section. The coupling 
conventions are discussed in detail in the text. The cross-sections are compared against the SM cross-section 
which yields 6.3 pb at the LO; the color scales for the models indicate the deviation from the SM cross- 
section, as indicated by the legend at the bottom. The curves indicate constant x^ for a given cross-section, 
as defined in Eq. (|22]). Contours for four cross-section values (cyan, blue, green, purple) are shown for 
X^/d.o.f. = 1 and 2. A single (cyan) x^/d-o.f. = 3 contour with SM cross-section is shown. Model points 
of a given color should be compared to x^ contours of the same color. 
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Detector 


CDF 


ATLAS 


{rj, 0) cells in cal 


(80,24) 


(81,63) 


r) width of cal cells for r? < 5 


0.1 


0.1 


(f> width of cal cells 
EM cal resolution (GeV) 
had cal resolution (GeV) 


0.262 


0.1 


0.01 e 0.2v^^/GeV 


0.01© 0.1 V^/GeV 


0.8y/E/GeV 


0.8 ^^/GeV 


MET resolution 


0.2 


0.2 


cal trigger cluster threshold 


3 GeV 


3 GeV 


outer radius of tracker 


1.0 m 


1.0 m 


magnetic field 


1.4 T 


2T 


sagitta resolution 


4 X 10"^ m 


5 X 10-6 m 


track finding efficiency 


0.98 


0.98 


minimum track Pt 


0.30 GeV/c 


0.3 GeV/c 


tracking eta coverage 


2.0 


2.5 


e/7 eta coverage 


2.0 


3.0 


fi eta coverage 


2.0 


2.4 


r eta coverage 


2.0 


2.0 



TABLE VL Detector Parameters of PGS4 simulation for the Tevatron CDF and LHC ATLAS detectors 

employed the cone algorithm, the LHC analyses use the anti-fcr-jet algorithm for jet reconstruction 
UTTII . We modified the kT-jet algorithm implementation in PGS by changing particle -particle and 
particle-beam distance measures according to the anti-fc^-jet algorithm definition. For comparison 
with Tevatron results, we use the cone algorithm with AR = 0.4, and for comparison with LHC 
results, we use the anti-fc-r-jet algorithm with AR = 0.4. 

After parton showering/hadronization and detector simulation which includes calorimeter er- 
rors in the measurement, the jet energy obtained from the jet reconstruction algorithm will be 
significantly different from the true value of the original parton. Note that this can be a significant 
source of distortion in event distributions with respect to energy scale variables such as invari- 
ant masses or transverse momenta. Therefore, such "measured" values of the jet energy must 
be corrected by performing standard candle experiments. The experiments have published their 
jet energy scale correction procedure in the literature, and we carry out our own jet energy scale 
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FIG. 23: Jet energy scale pseudo-experiment using e+e -^ qq process. 

correction for the PGS detector simulation. For this purpose, we generated the SM dijet event 
samples from electron-position coUisons (e+e~ — t- qq) with y^ = 200 GeV to 3000 GeV in 200 
GeV increments for the PGS implementation of the ATLAS detector, with 100,000 events for each 
energy sample. We generated similar samples for the PGS implementation of the CDF detector 
from v/s = 200 GeV to 1000 GeV in 100 GeV increments. 

From the samples, we compared the mean value of measured jet energy scales to the designated 



parton energy, and extracted the variance of the probabilistic jet energy measurement. In Fig. (23 1, 
we present the jet energy shift in this pseudo-experiment. The resultant jet energy scale distortion 
is reasonably matched with that of real detectors. The jet energy scale correction for CDF is given 
by 

Apr 1-63 



Ar/ = 0, 



(Bl) 



PTobs V^PTobs 

where probs is the nominal "p-r value of a constructed jet and Ap^ = PTtrue — PTobs is in GeV. For 
ATLAS, 

Apr AE 14.23 + 7.53r/2^3 Am 



E. 



Ari = 0, 



(B2) 



^^°^^ """^^ V^^obscosh^r/obs + m^,^ ^°b« 

where mobs is a jet mass. Here, for the CDF analysis, we ignore the jet mass, and jet momentum is 
parameterized only by p^, r], 0, while we retain a non-zero jet mass for the LHC analysis, since 
jet mass is a variable used in mu reconstruction. The variance of jet energy and angular parameters 
from SM dijet simulation for CDF is: 

a{pT) 1.21 _ 0.65 



0.0593 



PTohs 



y/PTohs ' 



a{r]) = 0.0112 



i/PTobs 



(B3) 
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FIG. 24: Cross-section versus x^ of top pair reconstruction for simulated LO SM events at ^/s = 1.96 TeV. 
The bin size below x^ = 30 is in increments of 3, corresponding to the number of degrees of freedom in 
the x^ fit. 

b. Top Reconstruction: To extract ApB, we reconstruct particle momenta of a semileptonic 
top pair using the x^ method, as in the CDF analysis. We find the missing neutrino momentum 
and fix combinatorics by minimizing x^ of over-constrained on-shell mass relations: 

yi=pl = 0, y2 = {Pe + Puf - m'^ = 0, ^3 = (pb^ +Pi+ Pyf - ml = 0, (B4) 

1/4 = {Pji + Pjif - m^ = 0, 2/5 = {Ph + Pj2 + Phuf - m^ = 0. (B5) 

X^ is defined hy x^ = y^ ■ V~^ ■ y where y = {yi,y2, 2/3, 2/4, 2/5) and V is the covariance matrix of 
2/jS. The detailed method is presented in Appendix A of [[6TII . 

Since we use a modified jet energy scale for the analysis, jet energies and the covariance matrix 
V are corrected correspondingly. V depends on individual jet-jet covariance matrix, which is 



related to a{pT) and a{ri) obtained in Eq. (B3 1. 



As a consistency check of our jet energy scale correction and top reconstruction algorithm, we 
show the resultant x^ distribution of SM tt events in Fig. (24). Taking into account degradation 



due to the combinatoric background, the result is reasonably well-matched with a theoretical curve 
and with the CDF analysis (shown as Fig. 15 in [|7]|). 

2. Event Generation and Parameter Scan Strategy 

We analyzed eight classes of models : CIV, CSV, F8C1V, CIS, CSS, schanC8Vr, schanC8VA, 
schanCSVR, as discussed in the main body of text. We consider the ti pair production cross 
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Model 


Parameter Scan Range (mass in GeV unit) 


CIV 


Rough: {im,gn)\m£ {200,400,600,800, 1000}, 5^ G {0.5, 1.0 ... 5.0} 

Fine: (m, gn) = (200, {0.4, 0.45 . . . 0.95}), (300, {0.4, 0.45 . . . 1.30}), 

(400, {0.6, 0.65 . . . 1.40}), (600, {1.0, 1.05 . . . 1.90}), (800, {1.30, 1.35 . . . 2.2} 


CSV 


Rough: {{m,gii)\me {200,400,600, 800, 1000}, g^ G {0.5, 1.0 ... 5.0} 

Fine: (m, gn) = (200, {0.2, 0.25 . . . 0.40}), (300, {0.3, 0.35 . . . 0.80}), 

(400, {0.4, 0.45 . . . 0.90}), (600, {0.5, 0.55 . . . 1.50}), (800, {0.7, 0.75 . . . 2.0}) 


F8C1V 


Rough: {(m, g, r])\m £ {200, 400 .. . 800}, 5 = 0.5, t] G {0, 0.5 .. . 3.0}} 
Fine: {(m, g, rj)\m G {300, 350 .. . 700}, 5 = 0.5, t] G {0, 0.5 .. . 3.0}} 


CIS 


Rough: {(m, gif)|m G {200, 400, 600, 800, 1000}, 5 G {0.5, 1.0 .. . 5.0}} 

Fine: (m, gn) = (200, {1.5, 1.55 . . . 2.0}), (300, {1.5, 1.55 . . . 2.20), 

(400, {1.5, 1.55 . . . 2.30}), (600, {2.0, 2.05 . . . 3.0}), (800, {2.5, 2.55 . . . 4.0}) 


C3S 


Fine: (m, gn) = (400, {1.5, 1.55 . . . 3.50}), (600, {2.5, 2.55 . . . 4.50}), (800, {3.5, 3.55 . . . , 5.5}) 


schanCSVr 


Fine: {{m,gR,n^,m^)\m e {A20,U0},gR G {0.35,0.45 ... 0.65}, n<^ G {4,5,6,7},m^ = 100} 


schanCSVA 


Rough: {(m,5|,fir*j)|mG {1600, 1800... 2400}, ff|j = -0.3, g*jG {1,2... 5}} 


schanCSVR 


Fine: (m, 5^, 5^) = (700, -0.05, {2.0, 2.5 .. . 6.0}), (850, -0.08, {2.0, 2.5 .. . 8.0}), 
(1000, -0.15, 3), (1000, -0.125, 5), (1000, -0.1, 8), (1500, -0.4, 5.5), (1500, -0.3, 8) 



TABLE VII: Summary of model points scanned. For the s-channel model with large decay width, 
schanCSVr, we take an additional contribution to the width of the mediator into scalars (p which is 

r^/m « (52n2/167r)(l - 4ml/m^f/^ EH- 

section at the LHC and Tevatron, and Aps at the Tevatron as the test of different models. We 
generate events for the process tt + or 1 jets with MLM matching. The renormalization group 
and factorization scales are fixed to be 200 GeV, and the top quark mass is 172 GeV. We employ 
CTEQ6L parton distribution functions. 

Our analysis has been done in three steps: rough scan, fine scan and benchmark point analysis. 
For the rough and fine scans, we generate 100,000 MadGraph events for each point, passing them 
through the PYTHIA and PGS pipelines. At the stage of the fine scan, we are able to see which 
models look most promising as shown in Figs. ([8|)-([T0]). For each benchmark model point, we 



generate five million MadGraph events followed by PYTHIA and PGS. Note that the number of 
generated events is reduced by 20% - 40% due to the MLM matching procedure. We summarize 
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the model points in Table [VII 



3. Cluster Pipeline Setup 

Although not a physics problem, generating event sets for a large number of model points is 
an intensive computing task with many engineering issues. The difficulty arises particularly with 
a cluster computing setup because it poses a new paradigm for software design. We share our 
experience in addressing such issues and suggest a common infrastructure. 

We utilize a cluster server named F lux in the Center for Advanced Computing at the University 
of Michigan. Ideally, the event generation software MadGraph could handle the cluster server 
configuration seamlessly, but in practice it is not easily implementable. The major challenge is 
that the disk I/O speed of a shared file system is not fast enough for MadGraph event generation. 
MadGraph creates a large number of small files, and the delay in writing to the shared file system 
causes the program to crash. A solution is found by making use of local storage; in most modem 
cluster computers, each node is provided with its own (fast) local storage space. To utilize this, 
the desired workflow must be to send a job which installs MadGraph on the temporary local disk 
space in the cluster node, followed by event generation tasks, uploading generated files, and finally 
erasing the temporary files. 

Since parallel computing is, in some sense, highly nondeterministic (on account of network 
latency or cluster usage traffic), jobs routinely fail. Therefore, it is important to make a highly 
resilient system for reducing the burden of bookkeeping of failed jobs. We design each job as 
a smart agent program which autonomously tests and monitors its own progression status. This 
requires us to make a central server for controlling assignment and checking the status of each job 
by having each job client report its status and wait for a new assignment for the next job if failed 
or finished. 

Many high energy physics programs contain legacy codes, and MadGraph is no exception. 
Due to incompatibilities or missing features, these codes often must be modified by users. How- 
ever, quick-and-dirty code repairs usually increase the complexity of a system. To control this, 
we make wrapping modules for external programs which is under our version control. By making 
those modules easily installable, the overall development becomes much simpler and easier in er- 
ror control. We call this system pipeline which is essentially a set of installable high energy 
physics program modules. 
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Job Queue Server 




Client side_(cluste_r node) 
Modules 





FeynRiiles 




MadGraph 




PYTHIA 








PGS 








Malhemalica 







storage 
Server 



Negotiation 




(a) Server configuration (b) Detail of client configuration and interaction 

FIG. 25: Cluster job queue system setup and interaction. 

For inter-process communication, we choose a standard web service interface, since HTTP 
protocols are not blocked in the usual firewall setup of a cluster. By standardizing the job speci- 
fication interface and each computer configuration, one can achieve flexibility and extensibility in 
routine high energy physics jobs. The job queue server retains information for each task for future 
documentation, and also effectively dispatches jobs. The web service choice has been superior in 
making a good user interface and utilizing common available tools. 

Fig. ( [25| ) shows our pipeline setup. We develop the system in haskell using Glasgow Haskell 
Compiler (ghc) 7.0. The job queue server is supposed to be always on and waiting for new jobs 
or new job requests from the client, which can run either in a cluster or on common desktops. If it 
runs on a cluster, a bootstrap script called clust eregg automatically installs a ready-made setup 
for a job client with MadGraph and the rest of the needed software. Since each job client sends 
its configuration when it requests a new job, the job queue server dispatches a new job for which 
the client is adequate (for example, if the job client does not have Mathematica, then Mathematica 
jobs are not assigned). A client also rechecks whether a job is doable with its current setup, and 
finally both parties handshake on the job assignment. After the negotiation, the job chent proceeds 
with the job according to the job specification from the server, and the job specification and high 
energy physics tools are interfaced with pipeline. After the job is finished, a job client sends 
its results to the storage server and wipes out the temporary files. Every step of the job status 
is reported to the job queue server for monitoring purposes. We will announce details of the 
pipeline and jobqueue systems elsewhere soon. 
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High performance computing facilities are now practically mandatory in high energy theory 
projects even for understanding the implication of the current state-of-the-art high energy experi- 
ments, especially in the LHC era. Building a common computing software infrastructure adjusted 
to high energy physics will harness our physics community in a very positive way. 
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